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ABSTRACT 

We report new observations of circumgalactic gas in the halos of early type galaxies obtained by 
the COS-Halos Survey with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. 
We find that detections of H I surrounding early type galaxies are typically as common and strong as 
around star-forming galaxies, implying that the total mass of circumgalactic material is comparable in 
the two populations. For early type galaxies, the covering fraction for H I absorption above 10 16 cm -2 
is ~ 40 — 50% within ~ 150 kpc. Line widths and kinematics of the detected material show it to be cold 
(T < 10 5 K) in comparison to the virial temperature of the host halos. The implied masses of cool, 
photoionized CGM baryons may be up to 10 9 — 10 11 M . Contrary to some theoretical expectations, 
strong halo H I absorbers do not disappear as part of the quenching of star-formation. Even passive 
galaxies retain significant reservoirs of halo baryons which could replenish the interstellar gas reservoir 
and eventually form stars. This halo gas may feed the diffuse and molecular gas that is frequently 
observed inside ETGs. 

Subject headings: galaxies: halos — galaxies: formation — quasars: absorption lines — intergalactic 
medium 



1. GAS AND QUENCHING IN EARLY TYPE GALAXIES |Bland-Hawthorn et aL|2007[ ), are consequences of ETGs 

The dichotomy of color, star-formation rate, morphol- 
ogy, etc. between spirals and early-type galaxies (ETGs' 
was established nearly a century ago, and has been h i ghly 
refined over the past decade (e.g . Baldry et al.|2004 |Bell 



et al.||2004| |Faber et al.||2007[) . Despite long-standing 
empirical measurements, the physical processes that di- 
vide the two populations are still uncertain. The main 
challenge is to identify the mechanism(s) that quench 
star-formation (SF) and that inhibited it over the past 
~ 10 Gyr. Dramatically different scenarios have been 
proposed for how star formation is quenched. Ejective 
feedback is generally caused by the galax y itself, and in- 
clude s super winds that expel gas (e.g. Springel et al. 
2005| or more general hea ting and stripping of the gas 
Tonnesen & Bryan 20091, as well as mergers. Con- 



living in massive halos. 

The paradigm of 'red and dead' galaxies as quiescent, 
passively evolving stellar systems is being reconsidered 
in light of findings that a non-negligible fraction of them 



versely, preventative feedback, such as shock-heating of 



the g as prior to its accretion (e.g. Dekel & Birnboim 



2006 1 or the disruption of in-falling clouds an d their dif- 



fusion into the galaxy's gaseous halo (Putman et al. 2011 
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he correspondi ng SFRs reach up to a few 
([Crocker et al.|201l|). Nevertheless, most ETGs 



are quenched, and remain so. These observations raise a 
new series of challenges: What is the source of the inter- 
stellar gas? On what time scale(s) did it arrive? Given 
this reservoir of cool gas, why does the SFR remain low? 

2. THE COS-HALOS SURVEY: MOTIVATIONS AND DATA 

Gas entering ETGs from their large-scale surroundings 
must pass through their circumgalactic medium (CGM), 
where it is detectable with sensitive UV absorption-line 
tracers (limiting column density Ah i > 10 13 cm~ 2 ). 
Quasar absorption line studies have attempted to link H I 
absorption (generally Lya only) and nearby galaxies, but 
the results have been inconclusive, suffering from small 
sample si zes and the inabili ty to control for galaxy type 
and SFR. |Chen et al.| ( |2001[ ) showed that Lya absorption 
is common around ga laxies of a wide range of morpho- 
logical types, whereas Chen & Mulchaey (2009) showed 



that strong H I absorbers (Ah i > 10 14 cm ' z ) do not cor 



relate with ETGs (see their Figure 13). Wakker & Savage 
( 2009 1 found ubiquitous H I around > 0.1 L* galaxies at 



350Ttpcscales, but could not separate galaxies by type. 
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Figure 1. Hi absorption profiles and SDSS images for two example passive (top) and star-forming (bottom) galaxies. The strong Lya 
absorption profiles consist of multiple blended components, and are similar for both galaxies. Individual Voigt profile components are 
shown in green; composite profiles are in red. The weaker Lyman series lines in the middle panels constrain the component structure and 
intrinsic line widths, which in turn constrain the gas temperature (§(5j. 

of ETGs occurs at much lower incidence than in star- 
forming galaxies, such that the presen ce of O VI out to 
150 kpc is related to star formation (Tumlinson et al. 
We will present the full analysis of the H I ab- 



Most recently, Prochaska ct al. (2011) found that strong 
Lya absorption (W r > 300 mA) arises in galaxy halos — 
including less massive ETGs — within p ~ 300 kpc, but 
the majority of their pairs lie at large radii (only 3 pairs 
at p < 200 kpc). Careful selection of QSO sightlines is 
therefore crucial to study the dependence of CGM gas 
within 150 kpc, and spectroscopy of the host galaxies is 
important to measure the current SFR and galaxy metal- 
licities. 

We designed the COS-Halos survey to probe the in- 
ner CGM of L* galaxies at z ~ 0.2 and to test the 
prevailing theoretical picture of "hot" and "cold" accre 
tion, which might be related to the observed star for 



mation dichotomy in galaxies (Keres et al. 2005 Dekel 
& Birnboim 2006). COS-Halos has several distinct ad 



vantages over previous work for assessing the amount 
and properties of CGM surrounding ETGs: (1) it sys- 
tematically covers impact parameters — 150 kpc and 
stellar masses log A/* = 10 — 11 with galaxies from both 
color-magnitude sequences selected prior to knowledge 
of absorption, (2) the wavelength range of COS and the 
galaxy redshifts (z > 0.1) generally allow for measure- 
ments of the weaker H I Lyman series lines that give 
better column densities and line broadening parameters 
than Lya alone, and (3) the galaxies are well charac- 
terized by ground-based optical measurements of their 
stellar masses, SFRs, and metallicities. COS-Halos has 
previously found that highly ionized O VI in the halos 



2011b] 

sorption for our complete sample of 50 galaxies in a fu 
ture paper (Tumlinson et al. 2012, in preparation). In 
this Letter, we present results for 16 ETGs from COS- 
Halos, showing that their associated H I gas is plentiful 
(§(3|, bound (§|ij), cold (§J5|, and may represent a large 
mass of baryons m their CGM (§ [6]). 

The COS-Halos survey obtainea 39 spectra of QSOs 
at z < 1 with the Cosmic Origins Spectrograph (COS; 
Green et al.||2012 ) onboard the Hubble Space Telescope 
(HST) under program GO-11598 (PI Tumlinson). These 
data were reduced, co-added , and continuum-normalized 
with procedures out lined by Mciring et al. (20TT|, Tum- 
linson et aL] ( |2011a[ ), and |Thom et al.| ( |201ip . We mea- 
sured column densities for the absorption associated with 
targeted galaxies using both the apparent optical depth 
technique (AOD; [Sembach fc Sav age 1992 ) and iterative 
Voigt pro file fitting with the meas ured COS line spread 
function (Ghavamian et al. 20091. Our measurements 
for H I column density are weighted averages of unsatu- 
rated Lyman lines or profile fits to damped absorption. 
Where all available Lyman series lines are saturated we 
set a lower limit to Nn i with the highest-order detected 
Lyman line. The multi-component Voigt profile fitting 
yields N, 6, and v for components within each system 
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Figure 2. Lyct rest-frame equivalent width W r versus impact 
parameter p for our COS-Halos sample of SF (blue) and ETG (red) 
galaxies. Note that only 15 ETGs appear in this figure, since in 
1 case we do not have coverage of the Lya transition. Even with 
four low values seen in the ETGs, there is no statistically significant 
evidence of a dependence on galaxy type (red vs. blue). 

that are used in § [5| 

The medium-resolution Keck and Magellan spec- 
troscopy used to derive the properties of the galaxies in 
th e sample (SFRs, m etallicities, etc) are described fully 
by Werk et al. (2012). The resulting sample contains 16 
galaxies classified as passive and/or early type, based on 
their lack of star-formation (sSFR = SFR/M* < lCT 11 

We caution that these fields 
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have not been uniformly surveyed for all galaxies down to 
low luminosity limits to eliminate all other possible asso- 
ciations, but these observations empirically describe the 
distribution of gas around ETGs as a function of impact 
parameter, whatever its source. Figure [I] shows an exam- 
ple of an ETG and star-forming galaxy from our sample. 
The composite image from the Sloan Digital Sky Survey 
(SDSS) is shown, along with the Lya absorption owing 
to the galaxy. Both galaxies have similar H I column 
densities and absorption profiles. 

3. STRONG H I IS COMMON AROUND EARLY-TYPE 
GALAXIES 

Figure [2] plots the Lya rest-frame equivalent width W r , 
as a function of impact parameter p, and shows the re- 
markable similarity in H I strength between the ETG 
(red) and SF (blue) samples. A KS test on W r does 
not provide significant evidence against the null hypoth- 
esis that the two distributions are drawn from the same 
parent population (D — 0.32, P — 0.19). The four Lya 
non-detections drive the difference in the KS statistic, 
and offer a hint that the distributions may be slightly dif- 
ferent, but we cannot separate the distributions in a sta- 
tistical fashion; a larger sample of ETGs will be needed 
to show if this difference is real. 

We also calculated the "covering fraction" /hit of H I 
detections. Even if we conservatively assume that the 
absorbers have the Nn i given by the lower limits for 
saturated lines, we find that the covering fraction above 
logA^Hi = 16 is /hit = 0.21 for SF galaxies, and /hit = 



0.38 for the ETGs. This colu mn density has signific ance 
as the threshold above which Stewart et al. (2011) pre- 
dicted essentially no absorbers surrounding halos at the 
ETG mass scale, owing to the transition between hot 
and cold accretion. Our covering fraction result rules 
out /hit = 0.01 at >99.9% confidence. The correspond- 
ing covering fractions at logA^i > 15 are /hit = 0.64 
(SF) and /hit = 0.56 (ETGs), again consistent to within 
the statistical errors of 0.15-0.2. Thus we find that strong 
H I absorption is common around the ETGs, and that to 
the limits of our survey, it does not appear significantly 
less often or weaker than in the halos of star-forming 
galaxies. 

4. THE STRONG H I IS BOUND TO ITS HOST GALAXY 

The detected H I shows a strong relationship to the 
targeted ETGs. It is found at separations consistent 
with being well inside the virial radius (i? V i r ~ 300 kpc), 
and it is closely associated with the stars in velocity 
space. Figure [3] shows the full velocity range of the de- 
tected H I absorption with respect to the systemic ve- 
locity of the galaxy's stars (set to v = 0). The relative 
errors between absorption and galaxy redshift (i.e. be- 
tween the wavelength solutions of the optical and COS 
spectroscopy) is ~ 30kms . The filled symbols mark 
centroid velocities of the fitted components, while the 
range bars show the full width at zero absorption, a 
measure of the greatest kinematic extent of the detected 
H I. These velocities are plotted with respect to the in- 
ferred dark matter halo mass, which is derived from the 
photometrical ly-estimated s t ellar mass using the mean 
relation from jMoster et ah (20101. The dashed curves 
show the escape velocity from the halos as a function of 
mass, at distances of 50, 100, and 150 kpc (from out- 
side to inside) . Nearly all the detected H I is found at 
velocities below the escape velocity. This is true even 
if we increase the velocity ranges by \/3 to account for 
unconstrained projection effects (since relative Doppler 
shifts are a sightline-projected lower limit to the true 3D 
space motion). It is also possible that some of the de- 
tected ionized gas is in fact bound to satellite galaxies in 
the vicinity of the targeted galaxjpj if so, the kinemat- 
ics suggest that the satellites are themselves bound and 
would be counted among the CGM mass budget for the 
host galaxy. We conclude that the detected HI is pre- 
dominantly bound to its host galaxy, and not escaping 
at the time it is observed. 

5. THE STRONG H I IS COLD COMPARED WITH HALO 

VIRIAL TEMPERATURES 

Because COS-Halos typically covers multiple Lyman 
series lines at z > 0.2 (Figure [I]), we can often derive 
line-broadening parameters that are not available when 
relying solely on Lya (which is almost always saturated). 
We perform profile fits to derive iVjj i and Doppler b for 
each component, where possible. These line widths pro- 
vide robust upper limits on the gas temperature under 
the assumption that the broadening is purely thermal, 
b = \p2MTjmZ. From the measured Doppler parame- 
ters (Figure [4]) we can firmly establish temperature up- 
per limits of T < 2 x 10 5 K (6 t herm < 60 kms -1 ), with 
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Figure 3. The full velocity range occupied by the detected HI 
absorption, with respect to the galaxy systemic redshift (v = 0), 
plotted versus the inferred dark-matter halo mass for our sample. 
The fitted component positions are indicated by the diamonds, 
while the full velocity ranges are shown by the vertical red bar. 
The three dashed lines mark the halo escape velocity at 50, 100, 
and 150 kpc, from outside to inside. All the H I detected in the 
ETGs is well inside the halo escape velocity, suggesting that is 
bound to the host galaxies. 



80% at T < 10 5 K (6 thcrm = 40 kins" 1 ). These line 
widths are upper limits based on the observed profiles, 
and can only decrease if the observed profile turns out to 
be composed of narrower unresolved components, or is 
significantly affected by no n-thermal broadening ( which 
has shown to be common; |Thom fc Chen|[2008| |Tripp| 
et al.||2008p | 

We therefore conclude that the detected H I traces 



mainly "cool" (possibly photoionized) gas, at tempera- 
tures well below the > 10 6 K virial temperatures of halos 
in this mass range. These temperatures resemble the ex- 
pectations for the "cold mode" of galaxy gas accretion. 
Together with the finding that the detected gas is plenti- 
ful and bound, this finding suggests that the deposition 
of cold gas into galaxy halos (from the inside or the out- 
side) and/or the formation of clouds at -C T v - lr are not 
fully quenched even in galaxies that have ceased to form 
stars. 

6. IMPLICATIONS FOR THE CGM MASS 

The total mass of cool gas in ETG halos is of interest 
in comparison to their stellar masses and interstellar gas 
reservoirs. Having measured the covering fraction and 
typical column density, we can estimate the detected H I 
mass from physical arguments and simple models. The 
simplest possible mass estimate comes from multiplying 
a surface density by an area: 
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11 There may also be shallow "broad Lyman a" (BLA) absorp- 
tion hidden in the profiles; we cannot draw any conclusions about 
the presence or absence an additional hot gas component in ETGs. 



Figure 4. On a component-by-component basis, the TVjj j and 
Doppler b parameter for H I around ETGs. Even if we assume 
that there are no non-thermal contributions to the line-broadening, 
most of the detected H I has implied temperatures of T < 2 X 10 
K. 

where we have taken a typical Nm = 10 16 cm -2 , at which 
the hit rate is /hit = 0.5 for ETGs, and p = 150 kpc. 
This mass of H I is strictly a lower limit because we have 
taken the typical minimum iVn i permitted by the pro- 
files of saturated lines and surveyed a radius that may 
not e ncompass the full CGM mass (c.f. Prochaska et al. 
20111. The corresponding total mass of CGM hydrogen 
is JWh = Mhi//hi, where /hi is the typical neutral frac- 
tion in the gas. Thus the ionization correction is the 
major factor in setting the total mass, but the ionization 
conditions, and thus the total mass of the CGM gas, de- 
pend on the unknown temperature and density of the 
detected material. Modestly overdense gas in the CGM 
should be exposed to the cosmic ionizing background ra- 
diation, which will further suppress the neutral faction 
and imply a much larger mass. 

To estimate what CGM masses are possible for the 
ETGs we show a simple halo model in Figure [5] that uses 
a power-law gas density profile exposed to a uniform ion- 
izing background and held at fixed temperature. The 
gas density follows a power-law dependence with radius, 
n oc (i?/i? V i r )~ Q , with a = 2, normalized to a specific cos- 
mic overdensity <5o at R V i r . The photoioni zation equilib- 
rium solutions were derived using Cloudy ( Ferland et al. 
19981 models and including ionization by the Haardt &; 
Madau photoionizing background (for z = 0.2) and the 
temperature-dependent collisional ionization equilibrium 
(CIE). We neglect any photoionizing contribution from 
the nearby ETGs since their star-formation rates are not 
high enough to contri bute significantly to the 1 Ryd back- 
ground in their halos ( Tumlinson et al.|2011a 
sities of 5q = 20 — 50 are required to match t 



Overden- 
le ETG H I 



detections with logTVui ~ 16. This profile gives a total 
mass of logAfn ~ H- The lower model adopts So = 2, 
and is a better match to the non-detections. Note that 
the results depend less on the chosen temperature than 
on the density normalization. From these simple models 
we conclude that the detected CGM is consistent with 
a mass Mqgm ^ 10 9-11 M Q . This mass is significant 
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Figure 5. TVjj j versus impact parameter p for all 16 ETGs. Lower 
limits are systems for which all available Lyman lines are saturated. 
The model profiles show two simple halo models for overdensity 
normalization 5 v i r = 50 (upper pair, green) and 5 v i r = 2 (lower 
pair, purple). Within each set of curves the lower curve corresponds 
to T = 10 5 K, while the upper bounding curves are for T = 10 4 ' 3 K; 
these are the two relevant temperature limits from the kinematics 
analysis above. The total gas masses inside 150 kpc are log M (<150 
kpc)= 9.6 - 11. 

by comparison with the interstellar gas masses detected 
in other ETGs (§1). At the high end these masses are 
potentially significant reservoirs of "missing baryons" in 
the budgets of these massive halos. We emphasize that 
these simple models are not rigorous derivations of the 
ionization correction, from which the actual gas column 
densities Nh can be inferred. We will publish separately 
a full analysis of the gas column densities based on mod- 
els calibrated to the many available metal lines in these 
absorbers (J. K. Werk et al. in preparation). Rather, 
these simple estimates show that ETGs harbor a signifi- 
cant mass of cold, bound gas in their CGM. 

7. IMPLICATIONS 

The COS-Halos survey has detected significant quan- 
tities of cool, likely photoionized gas bound to the ha- 
los of massive early-type galaxies. The cool gas does 
not appear to be correlated with star-formation, in stark 
contrast with O VI, which is strongly de pendent on 
the presence or absence of star-formation (Tumlinson 
et al. 2011b). The cool gas reservoir may contribute 



> ifj y of gas to the halos of ETGs, comparable to 

the interstellar reservoirs (§ 1), and exists at well below 
the halo virial temperatures. This material is typically 
well within the escape velocities of the host galaxies, so 
it may be re-accreted by the galaxies and eventually used 
for star- format ion. The question then becomes, why does 
the gas not do so? 

These ETGs have been quenched, in the usual sense 
that they are no longer forming stars. These strong H I 
absorbers and the cool CGM they trace were predicted 
to disappear at the end of cold accretion, which sho uld 
not operate in halos at this mass ( |Stewart et al.|2011[ ). If 
Mcgm ~ 1O 1O M is typical of the star forming galaxies 
in COS-Halos as well as the ETGs, then the ETGs prob- 
ably have less mass in their CGM relative to their stellar 



and halo masses than the SF galaxies (at least in cool 
gas). This might be the signature of quenching that can 
suppress but not completely remove cold gas from the 
halo, or that the gas can reaccrete long after quenching. 
It is also possible that the quenching mechanism did re- 
move all the CGM gas, after which cold clouds re-accrete 
from the IGM or reform from thermal instabilities in the 
hot corona. 

Models that attempt to explain the formation and evo- 
lution of ETGs will need to account for this significant 
budget of cold gas in their halos. While it is possible 
that some of this gas is the direct product of ejection 
by supernovae while star-formation was occurring, it is 
perhaps more likely to be gas cooled from their hot halos 
or falling in from the IGM. In this case, preventive feed- 
back may act to slow its accretion, allowing only small 
budgets of interstellar gas to enter the ISM and form 
stars. Some of it (especially the stronger systems) may 
be bound to satellites that are themselves bound to the 
host, and thus present or future contributors to the dif- 
fuse CGM. In any case, our findings indicate that the 
quiescence of star formation in ETGs generally cannot 
be attributed exclusively to the quenching of gas in their 
circumgalactic medium. 
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